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Abstract. We present photoionization models of the 
high excitation gas in the Extended Narrow Line Region 
(ENLR) of NGC 1068. The ENLR line fluxes have been 
calculated allowing for attenuation of the central-source 
ionizing continuum as a function of distance from the cen- 
tre. Diffuse continuum emission from low density ENLR 
gas is included as an important secondary source of ioniza- 
tion. The observed high excitation emission further than 
25 arcsec from the centre of NGC 1068 can be fitted by 
photoionization models using a central-source luminosity 
of 3.6 X 10^4 erg s'^ between IQi^-S and IQi^^ Hz, with 
the continuum shape attenuated by nuclear gas with an 
integrated column density of N^ — 10^^ cm~^. The re- 
flected soft X-ray continuum from the attenuating gas 
could be responsible for about 10% of the observed, re- 
solved circumnuclear soft X-ray continuum extending out 
to 15 arcsec from the centre (Wilson et al 1992). 

Key words: galaxies: active - galaxies: nuclei ~ galaxies: 
individual (NGC 1068) - galaxies: Seyfert - modelling: 
galaxies 



1. Introduction 

Spectropolarimetric observations of the Seyfert 2 galaxy 
NGC 1068 show broad, polarized Fell and Balmer emis- 
sion lines in the nuclear spectrum (Antonucci & Miller 
1985; Snijders et al. 1986; Miller & Goodrich 1990; Miller 
et al. 1991; Tran et al. 1992). These results suggest that 
this galaxy actually has a Seyfert 1 nucleus which is ob- 
scured from direct view but is seen in light scattered from, 
and hence polarized by, obscuring material which may be 
arranged in the form of a disk or torus surrounding the 
central nucleus. 

Inspection of the ENLR spectra shows the existence 
of ionized gas emitting high-excitation [Nev]A3425 and 
[Olll]A5007 lines at an approximate distance of between 
15" and 50" from the centre of NGC 1068 (Evans & Dopita 
1986; Bergeron et al. 1989). Lower ionization lines such as 
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[Oii]A3727 and [Nii]A6583 have intensity ratios relative 
to Ha) which arc similar to those seen in Hii regions, 
and do not reflect the degree of ionization which might be 
expected based on the detection of the [Nev]A3425 line 
(Evans & Dopita 1986). Furthermore the spectra of the 
ENLR show a low intensity of [Oiii]A5007 and HeiiA4686 
lines compared to the intensity of [Nev]A3425. 

Evans and Dopita (1986) modelled the ENLR spectra 
in NGC 1068 with two emission components: Hii regions 
and highly-ionized low density gas. They also suggested 
that the continuum seen by the high-excitation gas has a 
turn-on energy (simulating a photo-electric cut-off due to 
absorption by intervening material closer to the centre) 
varying over the range 20-60 Ryd. This decreases the in- 
tensity of the HeiiA4686 and [Oiii]A5007 hues compared 
to the intensity of [Nev]A3425. Bergeron et al. (1989) used 
optical emission-line ratios to infer that the ionizing con- 
tinuum seen by the ENLR gas may be two orders of mag- 
nitude larger than the observed soft X-ray continuum. 

Nazarova (1995) modelled the high-excitation ENLR 
gas in NGC 1068 with a central power-law continuum and 
an extended source of ionization. The extended source was 
assumed to be stars of temperature 80 000K located 1- 
2 kpc from the centre (or 14-28 arcsec adopting 1 arcsec = 
72pc (TuUy 1988)). The central region of NGC 1068 does 
have a ring of very luminous Hii regions 13 arcsec from 
the centre (Snijders et al. 1982; Bruhweiler et al. 1991). 
Adding an additional stellar source changes the shape of 
the incident spectrum which ionises the ENLR gas by 
adding a bump at 2.5 Ryd (corresponding to the peak 
emission energy of a star with temperature 80 OOOK) . This 
bump changes the population of the ionization stages. The 
0+^ moves to the higher stages 0+^ and 0+^ and He+ 
to He"''+. Similarly Ne+^ moves to Ne+^, although this 
change is small due to the higher ionization threshold. 
Overall, adding the stellar component leads to stronger 
emission in [Nev]A3425 compared to that of [Oiii]A5007 
andHeiiA4686. 

Ahhough modelling the ENLR spectra of NGC 1068 
with an additional stellar ionization source helps alleviate 
the "energy deficit" problem, there still exists a problem 
connected with the geometry of the extended region. Since 
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there are many galaxies with ENLRs which show strong 
emission in the high-excitation [Nev]A3425 Hne (Binette 
et al. 1996), it seems unlikely that they all have a strong 
contribution from a ring of very luminous H ii regions at 
just the right location around the nucleus. 

Recently, Binette et al. (1996) proposed that the high- 
excitation line strengths in ENLRs, together with other 
problems such as low electronic temperatures and the 
small range in the Heii/H/? ratio, could be successfully 
solved by assuming that the ENLR contains a combina- 
tion of matter-bounded and ionization-bounded clouds. 
These two cloud populations have different spectra and 
therefore can be combined to reproduce a wide range of 
observed ENLR spectra. A principal result was that the 
ionization-bounded clouds are photoionized by radiation 
from the central-source which has been attenuated by the 
matter-bounded component. 

In this paper we explore the advantages of such multi- 
component gas models in attempts to fit the observed 
ENLR spectrum of NGC 1068, adopting the shape of the 
central continuum as given by Pier et al. (1994). We also 
examine the effect of different central luminosities and dif- 
ferent amounts of attenuation of these continua on the 
predicted ENLR spectrum. Models of the observed ENLR 
emission are discussed in Section 2. We examine in Section 
3 if the model which best fits the ENLR line emission can 
also explain the observed extended soft X-ray continuum. 
The conclusions are given in Section 4. 



2. Modelling the ENLR spectra. 

2.1. The ionizing continuum 

The shape of the central-source continuum in NGC 1068 
has been derived from two sources. The UV and X-ray 
continua ware assumed from Pier et al. (1994), while the 
rest of the continuum shape was taken to follow that of the 
canonical AGN continuum given in Mathews and Fcrland 
(1987). These two continua arc shown in Fig.l. 

The bolometric flux from the scattering material in 
the spectral range lOi^^-^-lO^*^-^ Hz is 1.5 x 10"^° erg 
cm~^ s^^ (Pier et al. 1994), corresponding to a lumi- 
nosity of 3.6 X 10'*^ erg s^^ for an adopted distance to 
NGC 1068 of 14.4 Mpc (TuUy 1988). Adopting this as 
a minimum central-source luminosity, our photoioniza- 
tion models were calculated within the luminosity range 
3.6 X 10*2 to 3.6 X 1044 erg s~^ (Log L = 42.55-44.55). The 
luminosity upper limit was chosen to be approximately 
half the IR luminosity of NGC 1068 {Lir = 7.3 x lO^" erg 
s"-'^) found by Telesco and Harper (1980), as about half of 
the IR emission could be due to dust heated by stars. 



2. 2. Location of the [Nc v] zone 

Evans and Dopita (1986) and Bergeron et al. (1989) de- 
tected extended high-excitation emission from two knots 
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Fig. 1. The adopted continuum of NGC 1068 compared 
to the canonical AGN continuum of Mathews and Fer- 
land (1987). Both continua arc plotted in units of Photons 
Ryd^^ cm^^ s^^ with arbitrary offsets. 



in NGC 1068 located between 29.4 and 43.5 arcsec and be- 
tween 24.7 and 40.5 arcsec north-east of the centre. Apart 
from these, there is also a small [Nev] emission patch in 
the south-west between 18.2 and 28.8 arcsec from the cen- 
tre. However the strong [O ill] emission extends further out 
than the Balmer emission, and is detected up to 45 arcsec 
from the centre (Bergeron et al. 1989). In our calculations 
we have adopted a location for the high-excitation ENLR 
region in NGC 1068 (henceforth called the [Nev] zone) 
between 25 and 45 arcsec from the centre (or between 1.8 
and 3.24 kpc). 

2.3. Method 

All calculations were carried out using CLOUDY (version 
c90.01; Ferland 1991) for a plane-parallel slab of gas with 
solar abundances. Because the high-excitation emissivity 
in the ENLR depends strongly on the shape of the ion- 
izing continuum, we calculated the line fluxes taking into 
account the attenuation of the ionizing continuum by gas 
at radii between 5 and 25 arcsec from the centre (Sec- 
tion 2.4; Nazarova et al. 1997). The transmitted portion 
of the central-source continuum and the diffuse continuum 
(emitted by the low density gas) were calculated for 6 dif- 
ferent luminosities of the central-source within the lumi- 
nosity range given in Section 2.1, with 5 different column 
densities and different corresponding gas densities. The 
sum of the transmitted continuum plus the diffuse con- 
tinuum was then taken as the incident continuum at the 
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Fig. 3. The electron temperature profiles for a homoge- 
neous envelope of low density gas in the [Nev] zone for 
different luminosities of the central-source. 



starting radius of the [Ne v] zone (25 arcsec from the cen- 
tre). The shape of the unattenuated, central-source contin- 
uum and the sum of the transmitted and diffuse continua 
are shown in Fig. 2. The shape of the attenuated continua 
depend on both the luminosity of the central-source and 
the column density of the attenuating gas. The attenua- 
tion is large in the EUV/soft X-ray range when the gas 
is effectively ionization bounded for the dominant absorp- 
tion species. 

It should be noted that the shape of the sum of the 
transmitted and diffuse continua is quite different from the 
transmitted continuum alone. Adding the locally-emitted 
diffuse continua reduces the depth of the gap in the contin- 
uum between approximately 1 and 100 Ryd. The emission 
of the high-excitation [Nev]A3425 and [Oiii]AA4959, 5007 
lines depend on the shape of the continuum in this region 
as it controls the parent ion populations. The electron 
temperature of the [Nev] zone is shown in Fig. 3 for differ- 
ent luminosities of the central-source. The sharp drop in 
temperature after column densities 10^^ cm~^ and 10^^'^ 
cm~^ respectively for low luminosities (Log L = 42.55- 
42.85) is not shown. For high luminosities (Log L = 43.85- 
44.55) the temperature remains high enough to produce 
strong emission in the high-excitation lines even for high 
column densities (N/j = 10^2-23 ^ni^^). 

2.4- Attenuating gas geometry 

Since we use in our calculations the sum of the trans- 
mitted and diffuse continua as the incident continuum for 



the [Ne v] zone, we need to consider the geometry of the 
gas which attenuates the central-source continuum. As the 
physical size of the region containing the attenuating gas 
is fixed (5-25 arcsec from the centre), varying its total 
column density N^, between 10^° and 10^'* cm~^ implies a 
variation in gas density between 0.023 and 234 cm^'^, as- 
suming the gas is homogeneously distributed and that the 
filling-factor remains the same. The starting radius for the 
region containing the attenuating gas has been taken to 
be 5 arcsec, as this is believed to be the maximum size of 
the region containing the material obscuring the central- 
source from direct view (Pier et al. 1994). The contribu- 
tion from scattered central continuum to the illumination 
of the [Ne v] zone depends on the optically thickness for 
electron scattering (Tscat)- The scattered continuum will 
also be absorbed in the gas with corresponding optical 
thickness (Tabsorb)- As Tscat is 10"^*^ times less than Tabsorb 
we can neglect the contribution of the scattered continuum 
to the shape of the ionizing continuum for the [Ne v] zone. 



2.5. Low density gas in the ENLR of NGC 1068 

The high-excitation [Nev]A3425 line and part, possibly 
most, of the [Oiii]A5007 and HeIIA4686 fines beyond 25 
arcsec are produced by highly-ionized, low density gas in 
the ENLR. The emission from low density gas in the [Ne v] 
zone corresponding to various total column densities is 
shown in Fig. 4 for different luminosities of the central- 
source. The two parallel lines show the observed ratios 
in the line intensities for different parts of the ENLR of 
NGC 1068 derived from Evans and Dopita (1986) and 
Bergeron et al. (1989). From Fig. 4 it appears that the 
observed line fluxes could be fitted with a luminosity of 
the central-source higher than Log L = 43.85. 

Adding emission from high density gas, which could ex- 
plain the observed low ionization line emission (discussed 
below), will increase the emission in the H/3, [Oiii]A5007, 
and HeIIA4686 lines, but makes little contribution to the 
emission in [Nev]A3425. Therefore, we rule out luminosi- 
ties below Log L — 43.85 as they fail to explain the min- 
imum observed flux in [Ncv]A3425 even with a filling- 
factor of 1. For luminosities higher than Log L = 43.85 
the filling- factor in the [Ne v] zone might be less than 1 in 
order to fit the observed line intensities. 

We note that the predicted [Nev]A3425 flux is lower 
than [Oiii]A5007 for all luminosities, in agreement with 
the observed F(A3425)/F(A5007) line ratio which varies 
between 0.26 and 0.58 in different parts of the ENLR 
(Bergeron et al. 1989). As can be seen in Fig. 5, a small 
change in column density produces a quite dramatic 
change in the F(A3425)/F(A5007) ratio. This is due to 
the approximately 1.7 times higher photon energy needed 
to produce [Nev] than [Olll]. Hence the former is more 
sensitive to the column. The observed F(A3425)/F(A5007) 
ratio can be fitted for a small range of "critical" column 
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Fig. 2. The unattenuated central-source continuum (labelled incident) together with the sum of the transmitted and 
diffuse continua for 6 different luminosities of the central-source. Units are Photons Ryd~^ cm~^ s~^, with an arbitrary 
constant added to each continuum for display purposes. 



densities from 10^"''' cm ^ to 10^^ cm ^ for the luminosity 
range considered here. 



2.6. Two- component model of the [Nev] zone 

The presence of strong emission in low ionization lines 
such as H/3, [Oii]A3727, [Oi]A6300, [Nii]AA6548, 6583, 
and [Sli]AA6717, 6731 indicates emission from low filling- 



factor high density gas, similar to Hii regions, in the 
ENLR of NGC 1068. This emission is in addition to that 
from low density gas. The emission from high density gas 
was therefore also calculated for the 6 different luminosity 
models. Because the most probable column density ob- 
tained from analysis of the low density gas is N/j = 10^^ 
cm~^ (Fig. 5), the sum of transmitted and diffuse continua 
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Fig. 4. The theoreticallluxes of the H/?, HeIIA4686, [Nev]A3425 and [Oiii]A5007 lines (in units of ergcm"^ s"^ at the 
ENLR) from low density gas in the [Ne v] zone versus column density for different luminosities of the central-source. 



for this column density was used to calculate the contin- 
uum incident on the high density gas in the [Nev] zone. 

Composite photoionization models for three luminosi- 
ties of the central-source are given in Tabic 1. The ob- 
served and predicted line fluxes in the [Nev] zone from 
1 cm^ of the surface (the surface brightness scaled to the 
distance of the [Nov] zone) are given together with the 
relative contributions of the high (HDG) and low (LDG) 
density gas (HDG/LDG) to the H/3 flux, their densities 



(Log Ne(LDG) and Log Ne(HDG)) and the fiUing-factor 
f for the [Ne v] zone. The filling- factor was estimated by 
comparison of the observed and predicted H/3 flux. Table 1 
shows that most appropriate luminosity for the prediction 
of the observed line fluxes is Log L = 44.55. 

Binette et al. (1996) note that all of the NLR photoion- 
ization models which fit the strong emission lines predict 
an [Oiii]A4363/[Oiii]A5007 ratio < 0.01, which is smaller 
than the observed ratio ~ 0.015. The electron temperature 
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Table 1. The logarithms of the observed surface brightness of emission lines scaled to the distance of the [Ne v] zone 
and the modelled lines fluxes (in units of erg cm~^ s~^) for three values of the luminosity of the central-source. 



A 


Ion 




Observed 


Log L= 44.55 


Log L = 44.25 


Log L = 43.85 


3425 


[Nev] 




-3.59 - 


-4.09 


-4.18 


-4.30 


-4.42 


3727 


[On] 




-3.29 - 


-4.10 


-4.10 


-4.22 


-4.76 


3869 


[Neiii] 




-3.80 - 


-4.50 


-4.28 


-4.10 


-4.35 


4363 


[Oiii] 


< 


-4.30 - 


-5.45 


-4.88 


-4.97 


-5.12 


4686 


Hen 


< 


-4.34 - 


-5.21 


-4.97 


-5.00 


-5.10 


4861 


H/3 




-3.63 - 


-4.43 


-4.11 


-4.13 


-4.06 


4959 


[Oiii] 




-3.77 - 


-4.34 


-3.54 


-3.53 


-3.59 


5007 


[Oiii] 




-3.08 - 


-4.07 


-3.08 


-3.08 


-3.13 


6300 
6548 


[Oi] 
[Nil] 


< 


-4.60 - 
-3.95 - 


-5.30 

-4.82 


-4.11 
-4.42 


-4.23 
-4.50 


-3.33 

-4.73 


6563 


BI3 




-3.28 - 


-4.99 


-3.61 


-3.70 


-3.59 


6583 


[Nil] 




-3.45 - 


-4.17 


-3.95 


-3.70 


-4.25 


6717 


[Sii] 




-3.96 - 


-4.70 


-4.11 


-4.22 


-4.28 


6731 


[Sii] 




-4.18- 


-5.07 


-3.94 


-4.00 


-4.10 


HDG/LDG 

Log Ne(HDG) 

Log Ne(LDG) 

Logf 










1 
4 
0.37 
-1.50 


1 
4 
0.37 
-1.33 


1 
4 
0.37 
-1.00 



HDG/LDG - The relative contributions from High and Low density gas to the H/3 flux; 
Ne(LDG) - The density of Low density gas in units cm" 
Ne(HDG) - The density of High density gas in units cm" 
f - The filling-factor. 
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Fig. 5. The theoretical and observed 

[Nev]A3425/[Oiii]A5007 line ratio versus column density 
for different luminosities of the central-source. 



obtained from the observed [Oiii]A4363/[Oiii]A5007 ratio 



is often higher than the equilibrium temperature of models 
calculated with densities less than 10^ cm"^. In our mod- 
els this line ratio is 0.0157 for a luminosity Log L = 44.55, 
consistent with the typically observed ratio. 

If we suggest a higher luminosity than Log L = 44.55 
for the central-source it increases the emission in the high- 
excitation lines and simultaneously decreases the emission 
in [Oi]A6300. However it seems unreasonable to increase 
the proposed luminosity of the central-source more than 
the observed IR luminosity of Log Lm = 44.88 (Tclesco 
& Harper 1980). Therefore we believe that the best fit for 
the observed ENLR spectra in NGC 1068 corresponds to 
a central-source luminosity of Log L = 44.55. 

3. The reflected continuum 

Investigation of the NGC 1068 nuclear region in soft X- 
rays (Wilson et al.,1992) shows 3 X-ray sources; unre- 
solved emission associated with the Seyfert nucleus (55% 
of the total soft X-ray flux); resolved, circumnuclear emis- 
sion extending from the centre to a distance of 15 arcsec 
(23%); and large-scale emission up to 60 arcsec (22%), 
which has a morphology and spatial scale similar to that 
of the starburst disk. As electrons or dust in the ENLR 
could reflect some of the central-source continuum into 
our line of sight (Heckman et al. 1995; Cid Fernandes & 
Terlevich 1995; Tran 1995), it is interesting to see how 
much of the circumnuclear soft X-ray continuum could be 
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Fig. 6. The observed, resolved continuum and model re- 
flected soft X-ray continua within 15 arcsec. The continua 
are plotted in Photons Ryd^^ cm"^ s~^ on an arbitrary 
scale. 



nuclear emission scattered from the proposed attenuating 
gas with column density N;; — 10^^ cm~^. 

We used CLOUDY to calculate the reflected contin- 
uum for different luminosities. We assume the reflected 
continuum is the continuum emitted from the illuminated 
face of the cloud back in the direction of the central-source 
(i.e. into 27r sr). For a plane-parallel slab the proportion 
of the reflected continuum in the direction to the observer 
could vary from ~ 0.5 (when the ENLR is located in 
the plane of the sky) to zero (when we see the nucleus 
through the ENLR). The optical images imply that the 
galaxy disk, and therefore presumably the ENLR gas, in- 
clination is 40°±3° (Brinks et al. 1997). This implies that 
the fraction of the reflected continuum in the direction of 
the Earth is w 0.4. 

The intensity of the reflected continuum depends on 
both the luminosity of the central continuum and the col- 
umn density of the scattering medium. As the column den- 
sity grows with distance from the centre we expect that 
the contribution of the scattered light should also increase 
with distance. 

The reflected continua calculated by CLOUDY for 
the various models are presented in Fig. 6. In each case 
the reflected continuum comprises the scattered central- 
source continuum plus the diffuse continuum emitted by 
the gas. If the central-source luminosity in NGC 1068 is 
Log L — 44.55 then about 10% of the observed, soft X-ray 
emission from the circumnuclear component extending to 
a distance of about 15 arcsec could be due to reflected con- 



tinuum from a column density N;i — 10^^ cni^^. The ma- 
jority of the extended soft X-ray continuum in NGC 1068 
would be due to a different origin, possibly the 10^-10^ K 
gas seen in the ROSAT HRI soft X-ray image. This gas 
could result from a hot outflowing wind driven by the pu- 
tative central hard X-ray source, or through shocks driven 
by the radio jets (Wilson & Elvis 1997). 

4. Conclusions 

The extended high-excitation emission observed between 
25 and 45 arcsec from the centre in NGC 1068 (the [Nev] 
zone) has been modelled with a combination of low and 
high density gas illuminated by a central ionizing source. 
Before illuminating the high-excitation gas the continuum 
is attenuated by gas located within 25 arcsec. The best fit 
model has a central-source luminosity between lO^**'^ and 
2Q18.4 jj^ of 3.6 X 10^'* erg s^^, a low density gas compo- 
nent with Log Ne = 0.37 and a high density gas compo- 
nent with Log Ng = 4. From analysis of the fluxes of the 
high-exitation emission lines, we find that the attenuat- 
ing gas between the centre and the high-excitation region 
has an integrated column density of N^ = 10^^ cm^^. 
The reflected soft X-ray continuum from this attenuating 
gas could explain up to 10% of the observed, circumnu- 
clear soft X-ray emission identified by Wilson et al. (1992) 
extending from the centre out to a distance of about 15 
arcsec. 
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